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Dynamics of double membrane films is investigated in the long-wavelength limit 5/1 < 1 (5 is 
the wave vector and h is the thickness of the film) including the over damped squeezing mode. 
We demonstrate that thermal fluctuations essentially modify the character of the mode due to 
its nonlinear coupling to the transversal shear hydrodynamic mode. The renormalization can be 
analysed under condition g <C 1 (where g ~ T/k, T is the temperature and k is the bending 
module). The corresponding Green function acquires as a function of the frequency lu a cut along 
the imaginary semi-axis. At qh > ^fg the effective length of the cut is ~ Tq 3 jr\ (where 77 is the shear 
viscosity of the liquid) . At qh < ~Jg fluctuations lead to increasing the attenuation of the squeezing 
mode: it is larger than the 'bare' value by the factor \j-Jg. 

PACS numbers: 05.20, 82.65, 68.10, 82.70.-y 



I. INTRODUCTION 



The most distinctive property of amphiphilic molecules is their ability to spontaneously self-assemble into aggregates 
of various shapes. Typically the molecules spontaneously self-assemble into membranes which are bilayers of a thickness 
of the order of a molecular length. Different lyotropic structures constituted of the membranes have generated 
considerable current interest (see the books [|l]-[| and the reviews ) . Films composed of two bilaycr membranes 
sandwiching a thin layer of a liquid are widely spread in the lyotropic systems. They play also an essential role for 
various biological processes (one can note the so-called flickering phenomena in erythrocytes or red blood cells). In 
the paper we will examine dynamic properties of such double membrane films. 

The main peculiarity of a membrane is its negligible surface tension. Indeed, the membrane is immersed into a liquid 
and consequently its area can vary. Zero surface tension is the equilibrium condition with respect to the variations. In 
the situation shape fluctuations of the membrane are determined by the bending elasticity, the corresponding energy 
is 0| 

where the integral is taken over the membrane which is considered as a two-dimensional object, K\,Ri are its local 
curvature radii, and k is the bending rigidity module. Corrugations of the membrane induced by the thermal noise 
lead to loosing the orientation correlation of the membrane pieces at separations larger than the so-called persistent 
length £p |9| which can be estimated as 

£ p ~ aexp(27TK/T) , 

where T is the temperature and a is the thickness of the membrane. The shape fluctuations of the membrane lead to 
the logarithmic renormalization of the bending module k, examined first by Hclfrich |hJ and later by Forstc r EJ , the 
correct renormalization-group equation was derived by Peliti and Liebler |p^| , Kleinert fl3|| and Polyakov |14| . The 
explicit form of the one-loop RG equation is 

dn 3T 
d£ = ' 

Here £ = \n(r/a) and r is the characteristic scale. As follows from the equation the role of the dimensionless coupling 
constant is played by the quantity 
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Note that ln(£ p /a) ~ g _1 . For real membranes g ~ 10~ 2 — 10 -3 and consequently we can treat g as a small parameter. 
The smallness of g means that there exists a wide range of scales r < £ p where thermal fluctuations can be treated in 
the framework of the perturbation theory. 

Corrugations of the membranes in a double film can be decomposed into undulation (or bending) deformations and 
the squeezing ones. The bending deformations are characterised by the displacement u of the film as a whole from 
its equilibrium position and the squeezing deformation is characterised by variations of the film thickness h (which is 
the separation between the membares). We will bel ieve that in equilibrium the film lies along the X — Y plane. Then 
in the harmonic approximation it follows from (1.1) that the energy of the film is 



H= I dxdy k(V 2 u) 2 + -{V 2 hf 



(1.3) 



where both u a nd h are treated as functions of x and y and V is the two-dimensional gradient. 

At deriving (1.3) we neglected the interaction between the membranes. First, one should remember about steric 
interaction associated with a certain restriction of accessible co nfig urations for one membrane in the presence of the 
second membrane Bj. The explicit expression for the energy is |D| 



Hster = dxdy 



3n 2 T 2 
128nh 2 



(1.4) 



Due to the interaction (|l.4j ) two membranes can be treated as independent only on scales smaller than g x l 2 h. 
Therefore (O) is the main contribution to the energy if 



qh> y/g, 



(1.5) 



where q is the characteristic wave vector. Second, we should take into account the Van der Waals interaction. We 
assume that the same liquid is inside and outside the film. Then the Van der Waals energy is U& 



Hvdw = dxdy 



Ha 2 
2nh 4 ' 



(1.6) 



where H is the Hamaker constant. We can neglect the energy in comparison with (1.3) if 

H /o\2 



In the following we will believe that the thickness of the film is large enough for the inequality 

H /a\ 2 



9 2 > 



ti /ay 



to be satisfied. Then (1.5) is the only restriction enabling us to treat the energy (1.3) as the main contribution to the 
film energy. 



II. DYNAMICS 



We will examine the dynamics of the double membrane film in the long- wavelength limit qh <C 1 where q is the 



wave vector of the eigen modes of the film. Note that the inequality qh <C 1 is compatible with (1.5) since jC 1 In 
the limit g/i<l one should take into account the following variables describing the dynamics: the velocity of the film 
v, the displacement of the film u, the film thickness h and the densities of both membranes since they are conserved 
quantities. We will be interested mainly in the squeezing mode associated with the relaxation of the thickness h. 

To find dynamical characteristics of the film one should solve the conventional hydrodynamic equations in bulk 
supplemented by boundary conditions on both membranes. In the linear approximation the problem was solved by 
Brochard and Lennon jljj ], they found the dispersion law of the squeezing mode 

nh 3 q 6 

uj = -i — , (2.1) 

24?7 

wher e to is the frequency of the mode and r\ is the viscosity of the liquid surrounding the membranes. At deriving 



(2.1) one assumed that in equilibrium the film is flat. Note also the dispersion law |_7| 
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(2.2) 



of the bending mode also found in the linear approximation. Note that th e d ispersion law (2.1) is correct only if 
to neglect the direct i nter action of the membares that is at the condition (1.5) whereas the region of applicability 
of the dispersion law (2.2) does not depend on the interaction of the membranes since they move in- phase i n the 
bending mode. The clastic modes associated with variations of the membrane densities are harder than ( [2.l|]2.2| ) |fl8| . 
Therefore the only effect of the elastic degrees of freedom at examining the squeezing mode is the incompressibility 
condition 



V a v a 



0. 



(2.3) 



Here and below we believe that all variables characterising the film are functions of x and y and assume that Greek 
subscripts run over x and y. 



We will consider the renormalization of the dispersion law (2.1) of the squeezing mode due to fluctuational effects. 
For the purpose nonlinear dynamical equations of the film should be utilised. In the long-wavelength limit qh <C 1 
the equations can be derived phenomenologically. The reactive (non-dissipative) part of the equations can be found 
using the Poisson brackets method (see |l^] and also |l8|]) whereas the dissipative part of the equations is expressed 
via kinetic coefficients. One should know the expression for the energy H. of the system for writing both contributions. 
Actually the expression only for one Poisson bracket will be needed for us: 



{j a (xi,yi),h(x 2 ,y2)} = h(x 1 ,y 1 )V a [8(x 1 - x 2 )5(yi - y 2 )} 



(2.4) 



where j a is the two-dimensional momentum den sity of the film. The expression (^4|) is characteristic of any scalar 
conserved quantity of a film ijTsfl , the expression (2.4) is motivated by the fact that the two-dimensional mass density 
of the film is ph where p is the three-dimensional density of the liquid. Note that j a ss phv a since we believe that the 
membrane thickness a can be neglected in comparison with the film thickness h. 



The dynamic equation for the thickness h has the standard form following from (2.4) 



d t h + V a (v a h) = TV 



,SH 
5h 



(2.5) 



where dt = d/dt and T is the kinetic coefficient. The second power of the gradient appeared in (2.5) since the equation 
should support the conservation law of the liquid inside the film and therefore the right-hand side of t he e quation 
should be a full derivative at any H. 

gets 



In the linear approximation we can neglect the sweeping term in (2.5) due to 

1) one 



Substituting the harmonic expression (1.3) for the energy TC into (2.5) and comparing the result with 



2.1 



h 3 /(12 V ) 



(2.6) 



Note that T is inversely proportional to the shear viscosity coefficient. The point is that the dissipation described by 
r comes from viscous motion of the liquid surrounding the double membrane film which is excited hardly at large rj. 
The dynamic equation for j a has the following form plfl 



dtja - {W,ja} = J a 



(2.7) 



where J is the momentum flow from the bulk to the film. Just the term supplies the main dissipation of the film 
momentum and therefore we neglected the internal viscosity. In the linear approximation [ plj 



J a = -2riqv a 



(2.8) 



where q is the non-local operator which is reduced to multiplying by the absolute value of the wave vector q in the 
Fourier representation. The Poisson bracket {7i, j a } can be reduced to the divergence of the symmetric stress tensor 
for any energy Ti . Actually only the contribution associated with the Poisson bracket ( [2.4| ) and created by the 
harmonic energy Ql.3| ) is relevant for us. Then the equation fl2.7| ) is written as 



dtja + - W Q V 4 /i = -2-qqv a . 



(2.9) 



We will not present here dynamical equations for the variables j z , u and the densities of the membranes. The reason 
is that the equations for j z and u describing the bending mode decouple in the approximation needed for us from the 
equations (2.5,2.9). Actually the equations for j z and u are the same as for a single membrane, the corresponding 
nonlinear equations can be found in |l8| and also in |2^,^l|. As to the equations for the densities of the membranes 
they need a separate consideration which will be presented elsewhere. Th e on ly role of the degrees of freedom at 
analysing the squeezing mode is reduced to the incompressibility condition (2.3). 
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III. RENORMALIZATION OF SQUEEZING MODE 



As is seen from ( p.l[ ) in the long- wavelength limit the squeezing mode is very soft. This is the reason why one 
anticipates that fluctuational effects related to the mode are relevant. The effects are associated with nonlinear terms 
in dynamic equations and can be examined in terms of the diagrammatic technique of the type first developed by 
Wyld who studied velocity fluctuations in a turbulent fluid. In the work p3] the Wyld technique was generalised 
for a broad class of dynamical systems. A textbook description of the diagram technique can be found in the book by 
Ma |24fl . The diagram technique can be formulated in terms of path integrals as was first suggested by de Dominicis 
p5f and Janssen (26). In the framework of this approach apart from conventional dynamic variables one should 
introduce also auxiliary fields conjugated to the variables. Then dynamical correlation functions of the variables can 
be presented as functional integrals over both type of fields: conventional and auxiliary. The integrals are taken with 
the weight exp(£T), where X is an effective action which is constructed on the basis of nonlinear dynamic equations 
of the system. 

Being interested in the renormalization of the squeezing mode of the double membrane film we will take into account 
only the variables h and v a and the corresponding auxiliary conjugated fields p and fi a . We should also remember 
about the incompressibility condition (2.3) and impose the analogous constraint V Q ^ Q = on the field /i Q . Then, 
say, the correlation function of the film thickness h is written as 



(hihz) — J VhVvtrVpVfi^. exp (iX) h\h,2 , 



(3.1) 



where the subscript 'tr' implies that in the Fourier representation we should take only transve rse to the wave vector q 
components of the fields v and f i. The explicit expression for the effective action figuring in (3T) can be found using 
the dynamical equations (|j| 2^6 ) . It can be written as the sum of the reactive and the dissipative parts X = 2" rcac +2diss 
where 



X rcac = J dtcfr ^pd t h + pv a V a h + fi a d t j a - ^ a V A hV a h^ , 
X diss = J dtd 2 r j-ir K pV 6 /j + iTT{X7p) 2 + 2rjnq{v + il» j 



(3.2) 
(3.3) 



The detailed derivation of the effective action for the problem can be found in ]2(J,|2l| . 

Let us introduce the designations for the pair correlation functions. Taking into account only the transverse 
components of the fields v and fi we can write 



(h(t,r)p(Q,0)) = 
(« a (t, r)^(0,0)) 
(h(t,r)h(0,O)) = 

(u a (*,?>/? (0,0)) 



dui d 2 q 



(2tt)3 
dui d 2 q 



exp(—iujt + iqr)G(u>, q) , 



(2tt)3 
du> d 2 q 



exp(—iujt + iqr) 



Taf3 



(2tt) 3 
dw d 2 q 



exp(— iut + iqr)D(uj, q) , 



exp(— iujt + iqr) 



Saf 



9a 9/3 



G tI (u),q) , 



D iT (u,q) 



(3.4) 



(3.5) 



The correlation functions (pp) and (fi/j.) are equal to zero (what is the general property of the technique, see e.g. H). 
The functions D and D tr determine the pair correlation functions of the observable quantities and the functions G, 
Gtr are response functions. Therefore, say, the function G(u>) is analytical in the upper t o half p lane. 

It is possible to formulate the diagram technique for calculating correlation functions (^4,3.5). The harmonic part 
of the effective action X = X rcac + Xdi ss determines the 'bare' values of the response functions 



G (u,q) = -- 



1 



Gtr,o(w,g) = -- 



1 



u> + iTKq 6 /2 ' ' ' phu + 2ir/q 

The values of the 'bare' pair correlation functions satisfy the relations 

1 



ImG=^-D. 
AT 



ImG tI = —D tI 



(3.6) 



(3.7) 
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which are consequences of the fluctuation-dissipation theorem. Besides the harmonic part the effective action I 
contains terms of the third order which determine t he thir d-order vertices which figure o n di agrams representing 
the perturbation series for the correlation functions (3.4,3.5). On e can c heck the relations ( |3.7j ) order by order and 
consequently they are true for the 'dressed' correlation functions (3.4.3.5). Note that 



(3.8) 



The relation (3J3) is a consequence of (3/7), analyticit y of G(u>) in the upper half-plane and of t he asymptotic law 
G(oj) « — u)^ 1 which is correct at large ui. Note that (3.8) could be obtained directly from (1.3) since the integral 



over frequencies is just the simultaneous correlation function. 

The analysis of the diagrams shows that they contain infrared logarithms related to the lines representing the 
correlation function D ([3.5]). The lines produce factors 



TL 

(V a h(t,r)Vph(t,0)) = —6 af} , 
Zttk 



(3.9) 



wher e L = ln\hg 1 / 2 /r] 
(|3.9[) can be found from (3 



and r 1 is determined by characteristic external wave vector of the diagram. The expression 
if to remember about the condition (1.5). The presence of the logarithmic contributions 



implies that the main renormalization of a correlation function like G(u>, q) is produced by the degrees of freedom with 
the wave vectors much smaller than q. Therefore we should extract from the diagrammatic expressions for G(u>, q) 
only the contributions corresponding to the interaction with the degrees of freedom. 

The program can be realized directly on the language of the functional integral. Let us separate the variables h, 
p, v, fi into fast parts (with wave vectors larger than q), basic parts (with the wave vectors of the order of q) and 
slow parts (with wave vectors smaller than q). At calculating G(u>, q) we c an fo rget about the fast parts and keep the 
interaction of the basic part with the slow part. Then we obtain from fl3.2jj3.3D 

' (3.10) 



X = J dtd 2 r ^pd t h + pv a m a + p a d t j a — —p a \7 4 hm a — T—p\7 6 h + 2rjfiqv^ 



where h, p, v, fj, denote the basic parts of the fields, m a is the gradient of the slow part of h and dots designate 
irrelevant terms. The action (|3.10| ) is of the second order over h, p, v, /j, and consequently the integrals over the fields 
can be taken explicitly. Since m varies only weakly on the length q^ 1 we get 



G(u;,q) = -((phuj + 2tr 1 q)A- 1 ) m , 

G tI (uj,q) =-((u + i K rq 6 /2)A- 1 ) m , 

A = (phuj + 2irjq){uo + i K Tq 6 /2) - Kq 4 m 2 tI /2 , 



(3.11) 
(3.12) 
(3.13) 



where 



= 8. 



'a/3 



and the designation (. . .) m means averaging over statistics of m. At calculating ( [3.11 3.12 ) we substituted j — phv. 
Actually the terms with ph can be neglected and we omit the terms below. 

At averaging in (3.11,3.12) the statistics of m can be regarded to be Gaussian. The point is that only sim ulta neous 
correlation functions of m enter the expressions and the functions are described by the harmonic function (1.3). The 
pair correlation function of m is equal to ( |3.9[ ) . Therefore 

(mtr> = 2^ ' 



and we find from (3.11) 



+00 



(3.14) 



We see that G as a function of the frequency to have the cut along the imaginary semi-axis which starts from 
u) = —iTnq 6 /2 and goes up to — ioo. The effective length of the cut can be estimated as Tq 3 /rj what is the new 
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characteristic frequency related to nuctuations. Let us compare the frequency with the position of the pole in the 
bare expression 



Tq 3 /y 
Tnq 6 



{qhf 



(3.15) 



We conclude that fluctuation effects dominate in the region g 1 / 2 < qh < g 1 ^ 3 . Now, we can justify neglecting phus in 
comparison with r\q in the above expressions. At qh ~ 1 



and at qh ~ Jg 



phwj [r\q) ~ pn/ (r] 2 h) ~a/li<l, 



phaj/{r/q) ~ pn/{ri 2 h)g 2 < 1 



Performing Fourier transform of (3.14) over frequencies one gets 



G(t,g) 



TL 



q 6 t 



-1/2 



exp 



(3.16) 



The expression (3.16) is correct for positiv e tim e t, for negative times G(t) = due to the causality principle since 
G is the response function. We see from (3.16) that in the fluctuation region g 1 / 2 < qh < g 1 / 3 there appears an 
intermediate power asymptotics i -1 / 2 which at large times t is changed by the exponential decay. That means that 
the squeezing mode is described by a nonlocal in time dynamic equation. 

The above assertion is correct for the wave vectors q > yfg jh. In the limit qh <C J~g we return to the local equation 



(2.5) but with the renormalised kinetic coefficient T. The quantity can be found if to integrate the weight exp(iX) 
over the degrees of f reedom with the wave vectors q > Jg jh. The main effect is related to the sweeping term in the 
effective action (3.2). Due to the integration over the degrees of freedom with the wave vectors q > ^fg jh the term 
iTT(Vp) 2 in (3.3) for the long-wavelength degrees of freedom is renormalised and we found for the renormalised value 
of I': 



r = -L I dtd 2 r(v { l.r)h(l. Drill. OlMO.O);. 



(3.17) 



where averaging is performed over the degrees of freedom with the wave vectors q > ^fg jh. Using the renormalised 
expressions for the correlation functions we get the estimate T ~ g~ x l 2 Y ^s> T. 



IV. CONCLUSION 



We demonstrated that fluctuations essentially modify the character of the squeezing mode due to its nonlinear 
coupling with transversal shear hydrodynamic mode. Fluctuation effects lead to non locality of the equation for 
the mode, the corresponding Green function is ( |3.16 ). The new characteristic frequency of the mode related to 



fluctuations is uj ~ Tq 3 jr\ (q is the wave vector) which remarkably does not depend on bending elasticity. It i s 



important to distinguish the characteristic frequency from the attenuation of the membrane bending mode (2.2) 
having the same q 3 dependence on the wave vector. Let us stress that the strong fluctuation effects are observed 
only for dynamics. Static characteristics are not influenced by fluctuations because of the smallness of the coupling 



constant ( |l.2| ). That is the reason why only the harmonic part of the energy (1.3) is needed for us. 

Strong dynamic fluctuations of h occur for the wave vectors q > ^fg jh. For smaller wave vectors nuctuations of 
h are weak. Nevertheless even for the wave vectors there is a memo ry of the region of stro ng fluctuations which 
is the renormalised value of the kinetic coefficient T in the equation ( |2.5|) : The 'bare' value (^Jj) is substituted by 
r ~ g ~ x l 2 Y r. Note also that to analyse the dispersion law of the squeezing mo de i n the limit q <€. yfg jh star ting 



from (2J5) one should take into account besides the energy (L5) also the steric (|l.4| ) and the Van-der-Waals (1.6) 



contributions to the energy. As a result one finds 

'9tt 2 T 2 lOHa 2 



-iTq 



64nh 4 nh 6 



Let us discuss possibilities to check our predictions experimentally. The membranes can be studied by a variety 
of experimental techniques. Lately laser "tweezers" are becoming a useful tool for probing dynamical properties of 
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membranes. The technique enables one to obtain direct information about amplitudes and characteristic times of 
dynamical fluctuations for different objects constituted from membranes. For details see the monography |27j] and 
recent experiments p9|~|3l|] . We can mention as well force apparatus measurements pl| making possible to investigate 
dynamical response for two very thin lamellar systems confined between the walls and the classical light scattering 
experiments. Because of relaxation of membrane fluctuations, the scattered light has a broadened spectral distribution 
compared with the incident light. Despite the broadening is small, the modern technique of light beating (intensity 
fluctuation spectroscopy) is allowed to obtain information about eigen-modes of the system. 

The conclusions concerning the renormalization of the squeezing mode in our opinion are interesting both in its 
own right and as a new test of the membrane fluctuations. 
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